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An  Ambiguity  Function  Independent  of  Assumptions 
About  Bandwidth  and  Carrier  Frequency 

D.  A.  Swick 

Electronics  Branch 
Sound  Division 

Abstract:  With  the  use  of  a  mean-square  difference  criterio  1  to  distinguish  echoes  from 
two  targets  at  different  ranges  moving  with  different  velocities,  an  ambiguity  funtdon 
is  derived.  The  concept  of  a  modulated  carrier  is  avoided,  and  the  attual  Doppler  effect 
of  time  compression  or  expansion  is  used,  rather  than  the  more  usual  approximate  -i  of 
constant  frequency  shift.  Thus  this  function  can  be  applied  to  the  very-low-  .«  ,:cy 
broadband  signals  sometimes  employed  in  sonar  systems.  It  reduces  to  the  Win,  .it  ard 
ambiguity  function  in  the  care  of  two  targets  of  nearly  equal  velocity,  and.  in  general,  in 
a  narrow-band  approximation.  Some  properties  of  this  ambiguity  function  are  discussed. 


INTRODUCTION 

The  ambiguity  function  has  been  used  by  Woodward  (1)  and  others  (2-7)  to  study  range 
and  velocity  resolution  of  radar  and  sonar  targets.  In  most  of  this  work,  the  Doppler  effect 
has  been  approximated  by  a  frequency  shift,  which  is  valid  for  narrow-band  signals.  Cahiander 
(8)  has  used  statistical  detection  theory  to  study  the  wideband  signals  u:>ed  by  bats  for  echo- 
location;  Kelly  and  Wishner  (9)  have  extended  the  theory  to  include  the  effects  of  high  velocity 
and  accelerating  targets.  Their  ambiguity  function  reduces  to  that  of  Woodward  for  narrow¬ 
band  signals,  if  terms  of  the  order  of  cxpf±2i<uoO  are  neglected. 

By  means  of  a  simplified  model  for  the  echo  of  a  signal  from  a  moving  target,  a  mean-square 
criterion  is  used  to  derive  an  ambiguity  function.  This  function  can  be  used  to  study  the  reso¬ 
lution  of  echoes  from  two  such  targets.  !t  is  equivalent  to  the  “signal  function”  of  Cahiander 
(8).  although  the  derivation  and  application  are  quite  different,  and  to  the  ambiguity  function 
of  Kellv  and  Wishner  (9).  Here,  however,  neither  bandw  idth  nor  carrier  frequency  enter  into 
the  derivation.  In  the  case  of  sonar  signals,  the  velocity  of  the  target  v  need  ~ot  be  very  small 
relative  to  the  velocity  of  propagation  of  the  signal  c;  all  that  is  required  is  t..at  v  <  c.  By  suit¬ 
able  identification  of  parameters,  the  new  ambiguity  function  can  be  identified  with  W'oodward’s 
function  in  the  case  of  two  targets  of  nearly  equal  velocities,  and  in  the  narrow-band  approxima¬ 
tion  to  the  detection  of  echoes  from  a  single  target. 

THE  DOPPLER  EFFECT  MODEL 

Let  s(t)  be  a  real-valued  function,  square  integrable  on  (— »,  «),  representing  a  signal 
transmitted  with  a  propagation  velocity  c,  assumed  constant.  Let  r(f)  be  the  radial  distance 
and  v  be  the  velocity  (assumed  radial  and  constant)  of  a  reflecting  object.  The  signal  portion 
of  the  received  (reflected)  waveform  is  then 

x(t)  =  a  s[t  —  T{t)],  (1) 

where  a  is  a  constant  (possibly  complex  to  account  for  phase  change  on  reflection)  and  7(t) 
is  the  time  required  for  the  signal  to  reach  the  object  and  return.  Thus,  a  signal  received  at 
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time  t  was  transmitted  at  time  t  —  T(t),  and  jvas  incident  on  the  target  at  time  t  —  T(t)l2,  when 
'die  target  position  was  r<>  +  vt  —  vT(t){2.  Hence, 


cT(t)l 2  =  r0  4-  vt  —  vT(t)/2 


or 


and,  from  this, 
where 

is  the  “Doppler  stretch  factor"  and 

is  the  delay  of  the  signal  at  t  =  0. 


no  -  2<?+ . 

(2) 

x(t)  =  as(at  —  To), 

(3) 

_  c  —  V 

c  +  t> 

(4) 

2r„ 

To  c  +  v 

(5) 

THE  AMBIGUITY  FUNCTION 

If  there  is  another  target  with  radial  velocity  v'  at  a  radial  distance  r'(t)  —  rd  +  v't,  the 
signal  reflected  from  it  will  be 


x'(t)  =a’s(a't-To) 


(6) 


where  the  constants  are  defined  as  before.  We  wish  to  distinguish  the  waveforms  given  by  Eqs. 
(3)  and  (6).  Hence,  following  Woodward  (1),  we  would  like  their  mean-squared  difference 


c 


[x(t)-x'(t)]*dt 


(7) 


to  be  as  large  as  possible  for  all  values  of  the  parameters,  except,  of  course,  in  a  small  region 
near  x(t)  —  x:(i),  when  the  targets  are  in  fact  indistinguishable. 

Expanding  Eq.  (7)  we  get 

I  [x(t)-x'(t)]*dt  =  J  x*(t)dt+f  x'*(t)dt- 2[  x(t)x'(t)dt 
=  a*f  s'(aC-To)  dt  +  a'*  f  s' (a' t  -  T'0)  dt 


-2  aa-f 


s(at-To)  s(a't  —  T'c)  dt 


^a'a'  +  aa'J  _2aa,  C  j(at _  To)  j(a-x _  r)  dt  , 


(8) 


where  we  assume  the  signal  to  be  normalized  so  that 
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Since  the  first  term  on  the  right  of  Eq.  (8)  does  not  depend  on  the  signal  waveform,  and 
since  a  and  a'  m.«y  be  complex,  we  see  that  we  can  maximize  Eq.  (7)  by  choosing  a  signal  which 
minimizes  the  modulus  of  the  second  term  on  the  t  ight  in  Eq.  (8). 

Equivalently,  we  can  define  the  correlation  function 

0(r,y)=f  s(t)  s{yt  +  r)  dt ,  (9) 

where  either 

y  =  —  and  r  =  ^-  To  -  To  (10) 

or 

y  =—,  and  r  =  A  T'0  -  jT0,  (11) 

a  a  '  ' 

and  require  that  |0(r,  y)|  be  as  small  as  possible  except  in  the  vicinity  of 

0(0,  1)  —  J  s'-(t)dt  =  1.  (12) 

Equation  (12)  represents,  of  course,  the  correlation  of  signals  reflected  from  targets  at  the 
same  range  and  with  the  same  velocity. 

If  we  define  the  ambiguity  function  as 

4>(t,  y)  4  |0(t,  y)  j*,  (13) 

ih'm  the  “dislinguishability"  criterion  requires  that  4»(r,  y)  be  as  small  as  possible  except  near 
t/t(0, 1)  =  !.  The  functions  0{r,  y)  and,  hence,  t/»(r,  y)  are  independent  of  any  assumptions  about 
bandwidth  and  carrier  frequency. 

THE  RELATIONSHIP  TO  THE  WOODWARD  AMBIGUITY  FUNCTION 

lx-t  s(t)  be  the  Hilbert  transform  of  s(f ).  The  "pre-envelope"  of  s(t)  is  defined  by  Dugundji 
(10) to  be 


z(t)  =  s(t)  +  i  s(f).  (14) 

Then  s(t)  —  /?e{z(f)},  and  Eq.  (9)  becomes 


6(t.  y)  —  [  Re{z(t)}  Re{z(yt  +  t)}  dt 

^  -  x 

=  (1,4)  f  [z(t)  +  z*(0]  [z(y t  +  r)  +  z*(yt  +  t)}  dt 

J  —  x 

=  ( 1/2)  /?ejj  +  t)  +  z(t)  z(y/  +  t)  d/j 

=  (1/2)  Re{6,(r.  y)  +  0z(r,  y)}. 


(15) 
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where  the  asterisk  denotes  the  complex  conjugate 


and 


Mt,  y)  =  J 
0j(t,  y)  =  J 


z(t)  z*(yt  +  t)  dt. 


z(t)  z(yt  +  r)  dt . 


(16) 


(17) 


If,  following  Kelly  and  Wishner  (9),  we  were  to  express  s(f)  as  a  modulated  earlier, 


5(0  =i?«{z(*)}  =/?e{[2(t)e-<*”'] 


(18) 


then  $i(r,y),  Eq.  (17),  would  be  related  to  the  terms  of  order  exp(±2/a>0f)  and  could  be  ne¬ 
glected  in  a  narrow-band  approximation,  as  is  done  L»  them.  We  do  not  need  to  do  this,  however, 
for  in  fact  8i(t,  y)  vanishes  identically.  If 


Z(w) 


-f. 


z(f  )e~,“r  dt 


is  the  Fourier  transform  of  z(f),  then  we  can  write 

1 


(t,  y)  = 


4  7T1 


Ilf 


Z(w)  Z(v)eity  *  v" u  eimT  dt  dv  dto 


=  y-  j  Z(u>)  Z( — yw)e'“T  dco  =  0, 


09) 


since 


and 


Z(a>)  =  0,  m  <  0 
Z(  — yw)  =  0,  co  >  0,  for  y  >  0. 


(See,  for  example.  Refs.  10  and  1 1.) 

We  can  write  s(t)  as  in  Eq.  (18),  where  co »  is  now  an  arbitrary  parameter  which  may  be,  but 
does  not  have  to  be,  associated  with  a  “carrier."  If  we  iet /(f)  =z(f)  exp[— icoot],  Eq.  (16)  becomes 

fl»(r.  y)  =  j  /(f)  e'-7*(y t  +  T)e-"°(y'  *  r>  dt 
=  e  '“*T  f  f(t)f*(yt  +  t)  e~,d'  dt 

*  -  x 

-e-'-Vxfr.  <M,  (20) 

where 

1  )«»o,  (21) 

and 

f(t)f*(t  +  T)e-u,dt  (22) 


NRL  REPORT  6471 


5 


is  the  combined  time  and  frequency  correlation  function  of  Woodward  (1).  The  approximation 
in  Eq.  (20)  holds  in  the  case  of  two  targets  of  nearly  equal  velcxits.  It  requires,  of  course,  that 
fit)  he  slowly  varying  as  compared  to  exp  )  —  exp[— i(y  —  l)a>0/],  but  it  must  be  remem- 
betcd  that  coa  is  arbitrary. 

If  we  consider  the  echo  from  one  target  (t.e.,  let  a'  —  a'  —  1,  To  —  0),  and  if  we  identify  o>0 
with  a  single  transmitted  frequency,  then  </>  is  the  actual  Doppler  shift  of  this  single  frequency. 
It  is  a  constant,  but  no  assumptions  have  Iteen  made  about  the  constancy  of  Doppler  shift 
across  a  band  of  f  requencies. 

From  Eqs.  (13),  (15),  (19),  and  (20),  we  hase 

■Mr,  y)  =  ( 1/4)  |/?e{0i  (t,  y)}|*  *£  (l/4)|«,(r,  y)|*  =  (1/4)  <f>) ,  (23) 

where  <t>  is  given  by  Eq.  (21),  and 


iMt,  <f>)  =  lx(T«  4>)il 


(24) 


is  the  Woodward  ambiguity  function. 

SOME  PROPERTIES  OF  THIS  AMBIGUITY  FUNCTION 

By  analogy  with  'he  properties  of  the  Woodward  ambiguity  function,  as  discussed  by  Sieberi 
(12)  and  others,  some-  of  the  properties  of  the  new  function  can  be  s  tmmari?ed. 

If 

S(oj)  =  |  s(t)e~>m‘  dt 

is  the  Fourier  transform  of  s(t),  we  can  write  for  Eq.  (9) 

®(t,  y)  =  j J J  S(v)  S(w)  e<(*'  +  ymU  ei,,T  dto  dv  dt 

=  ^  f  S(w )  S(— yco)  elmT  d<o 

=  —  J  S(w)  S*(yw)  e'"r  dto. 


since  s(t)  is  real. 

By  consideration  of 


J  [s(0  —  Xs(y<  +  t)]*  dt  2s  0 


for  all  X,  and,  in  particular,  Eq.  (26)  2*  0  for  X  =  Vy  sgn  8(t ,  y),  it  follows  that 

1 


0(t,  y)  «  —p,  for  all  r,  y 
Vy 


(25) 


(26) 


(27) 


and 


0(t,  y)  ^  1 ,  for  all  r,  for  y  »  1 . 


(28) 
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Integration  of  0*{t,  y)  over  all  values  of  its  parameters  must  be  limited  to  at  most  0  <  y  <  *, 
as  pointed  out  by  Cahiander  (8),  who  shows,  as  do  Kelly  and  Wishner  (9),  that  the  usual  volume 
constraint  on  the  ambiguity  function  mu^t  be  replaced  by 

f  f  #*(t,  y)  dr  dy  —  f  —  |S(w)|*do».  (29) 

*0  "  -  x  ■'o  *** 

Here  |S(to)|*  must  vanish  sufficiently  rapidly  as  to  —*  0  so  that  the  integral  in  Eq.  (29)  exists 
(9).  Equation  (29)  is  not  independent  of  signal  waveform  as  is  the  equivalent  integration  of  the 
Woodward  ambiguity  function.  In  a  narrow-band  approximation,  however,  the  integral  reduces 
to  ir/too,  a  constant  (8,9). 

In  place  of  the  twofold  symmetry  axis  of  the  Woodward  function,  we  have  for  y  ^  0 

0(~y'  y)  =  J  *(*)  *(~  “")  dt  =  y  J  s(t)  s(yt  +  r)  dt 

=  y 0(r,y),  (30) 

which  shows  the  relationship  between  the  alternative  definitions  of  the  parameters  r  and  y, 
given  in  Eqs.  (10)  and  (11). 

In  the  case  y  —  1,  Eq.  (9)  becomes 

0(t,1)  =  J  s(t)  s(t  +  r)  dt  =  R(r), 

*'-x 

an  autocorrelation  function  of  s(t)  ;  and  Eq.  (16)  becomes  an  autocorrelation  function  of  the 
“pre-envelope"  z(t). 

The  Woodward  ambiguity  function  is  related  to  the  output  of  a  filter  matched  to  a  signal 
s(f),  in  response  to  a  signal  s(t  —  r )  exp(— irf>t),  i.e.,  an  echo  with  a  constant  time  delay  and  a 
constant  frequency  shift  (13).  Similarly,  the  last  term  on  the  right  of  Eq.  (8),  and  hence  Eq.  (9), 
is,  apart  from  a  constant  multiplier,  the  output  of  a  filter  matched  to  s(at  —  T0)  when  the  input  is 
s(a't  —  T0)  (9),  evaluated  at  the  time  the  output  peak  signal  power  occurs.  If  we  let  a  =  a'  =  1  and 
T'f,  =  0,  then  0(r,  y)  =  6{~T0,  a)  is  the  output  of  a  filter  matched  to  the  transmitted  signal,  when 
the  input  is  an  echo  from  a  moving  target.  Each  frequency,  to,  in  the  transmitted  signal  has  been 
shifted  by  an  amount  proportional  to  itself,  i.e.,  is  received  as  &></= aa»,  where  a  is  given  by  Eq.  (4). 
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